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Although symmetry methods and analysis are a necessary ingredient in every physi- 
cist's toolkit, rather less use has been made of combinatorial methods. One exception is in 
the realm of Statistical Physics, where the calculation of the partition function, for exam- 
ple, is essentially a combinatorial problem. In this talk we shall show that one approach 
is via the normal ordering of the second quantized operators appearing in the partition 
function. This in turn leads to a combinatorial graphical description, giving essentially 
Feynman-type graphs associated with the theory. We illustrate this methodology by the 
explicit calculation of two model examples, the free boson gas and a superfluid boson 
model. We show how the calculation of partition functions can be facilitated by knowl- 
edge of the combinatorics of the boson normal ordering problem; this naturally gives rise 
to the Bell numbers of combinatorics. The associated graphical representation of these 
numbers gives a perturbation expansion in terms of a sequence of graphs analogous to 
zero-dimensional Feynman diagrams. 
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1 Introduction 

According to the axiomatic approach, it is possible to reconstruct a quantum 
field theory from the vacuum expectation values. When we normally order a set 
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of operators so that all the annihilation operators are to the right, then in such a 
vacuum expectation value only the constant term remains.^) Similarly, when we 
normally order a string of boson operators a,a^ satisfying [a, a^] = 1 , then the 
vacuum expectation reduces to the constant term; and likewise we may recover the 
expectation in the coherent state \z) by replacing a by the c-number z. Therefore 
the ability to normally order a string of operators is a powerful calculational tool in 
physics. It was recognized early on that the procedure of normally ordering a string 
of bosons leads to classical combinatorial numbers, the Stirling and Bell numbers 
Apart from the generation of interesting combinatorial sequences, and their 
extensions [2] , the associated graph theory representations of these numbers have 
analogies to the Feynman diagrams of a zero-dimensional quantum field theory. 

In this note we shall present a standard graphical representation of the Bell 
numbers which we show is in some sense generic. That is, we use these graphs to 
provide a representation of a perturbation expansion for the partition function of 
a free boson gas, and show that the same sequence of graphs can be used to give a 
perturbation expansion for a general partition function. 

We illustrate this approach by applying it to a supcrfluid boson system. 

2 Combinatorial sequences, generating functions and graphs 
2.1 Some combinatorial sequences 

Sequences of numbers of combinatorial interest are common in mathematics; for 
example, some sequences depending on one parameter are 

1. {n} The integers. 

2. {n\} The number of ways of putting n different objects into n different con- 
tainers (leaving none empty). 

3. {2"} The number of maps from an n-element set to {0, 1}. 

4. {B{ti)} The Bell numbers; the number of ways of putting n different objects 
into n identical containers (some may be left empty). 

These are all familiar except perhaps the last, the Bell numbers. These have the 
values 

B{n) = 1, 2, 5, 15, 52, 203, . . . , n = 1,2, ... . 

and B{n) < nl. 

We may also define sequences which depend on, say, two parameters, for exam- 

n\ 

pie the 'choose' symbols "Ct = -t — r. Related to the Bell numbers are the 

^ ^ (n-fc)!fc! 

Stirling numbers of the second kind S{n,k), which are defined as the number of 

^) Note that this process of normally ordering, denoted by J\f, does not change the value of the 
operator. Thus Af{f{a, a^)) = /(a, at). 
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ways of putting n different objects into k identical containers, leaving none empty. 
From the definition we have 

n 

B{n) = J2S{n,k). (1) 

fe=i 



2.2 Normal order 

Although somewhat unfamiliar to physicists, the Bell and Stirling numbers are 
fundamental in quantum theory. This is because they arise naturally in the normal 
ordering problem. For canonical bosons [a, a^] = f the Stirling numbers of the 
second kind S{n, k) intervene through 

n 

(atfl)" = ^5(n,fc)(at)'=a'=. (2) 
fe=i 

The corresponding Bell numbers B(n) = '^('^i ^) ^'^s simply the expectation 

values 

B{n) = (z|(ata)"|z},=i (3) 

taken in the coherent state defined by 

a\z)^z\z) (4) 

for z = 1. In fact, for physicists, these equations may be taken as the definitions 
of the Stirling and Bell numbers. 



2.3 Generating functions 

Corresponding to a sequence of combinatorial numbers {a„} we may define the 
exponential generating function (egf) 

A(x)=^a„— . (5) 
«=o ' 

This provides a compact formula for the combinatorial sequence. For example for 
the sequences above 

1. a„ = n, A{x) — xe^, 

2. a„ =n!, A{x) = 1/(1 - x), 

3. a„ = 2", A{x) = e2==, 

4. On = B{n), A{x) = exp (e"" - 1). 
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Except for the last, which we shall prove in what follows, these exponential gen- 
erating functions are immediate. Note that we consider the expansion in a formal 
sense; we are not concerned here about the convergence of the series defined by the 
egf 's. We may similarly define egf 's corresponding to 2-parameter sequences, thus: 



n=0 ^ k=0 ^ ' 



(6) 



For the 2-parameter sequences above we obtain: 

1. a„^fe="Cfc, = exp((l + 

2. anM = S{n,k), A(a;, y) = cxp(y(e^ - 1)) . 

In the latter function, the coefficient of the n-th power of x is an n-degree polyno- 
mial in y, a Bell polynomial, Bn{y). 

Therefore we have the generating function for the Bell numbers and Bell poly- 
nomials 



ABix) = J2B{n) 



n=0 



exp(a;+- + - + ...) (7) 



and 



AB{x,y) = Y,[T.S{n,k)y^)- = 



n=0 ^ k=0 
n=0 



exp(j/(e--l)). (8) 



2.4 Graphs 



We now give a graphical representation of the Bell numbers. Consider labelled lines 
which emanate from a white dot, the origin, and finish on a black dot, the vertex. 
We shall allow only one lino from each white dot but impose no limit on the number 
of lines ending on a black dot. Clearly this simulates the definition of S{n, k) and 
B{n), with the white dots playing the role of the distinguishable objects, whence 
the linos arc labelled, and the black dots that of the indistinguishable containers. 
The identification of the graphs for 1,2 and 3 lines is given in the Figure 1. 

We have concentrated on the Bell number sequence and its associated graphs 
since, as wc shall show, there is a sense in which this sequence of graphs is generic. 
That is, we can represent any combinatorial sequence by the same sequence of 
graphs as in the Figure 1, with suitable vertex multipliers (denoted by the V terms 
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n=1: O • V, 

S(1,1)=1 B(1M 




S(2,2)=1 S(2,1)=1 B(2)-2 




S(3,3)=1 S(3,2)=3 S(3,1)=1 B(3)=5 

Fig. 1. Graphs for B{n), n = 1, 2, 3. 

in the same figure). One reason the graphical representation is useful is that there 
exist some rather powerful results which apply to graphs and their associated gen- 
erating functions. We give one such result now. 

2.5 Connected graph theorem 

This states that if C{x) — c{n)x" /nl is the exponential generating function 

of labelled connected graphs, viz. c{n) counts the number of connected graphs of 
order n, that is with n lines, then 

A{x) = exp(C(a;)) (9) 

is the exponential generating function for all graphs 31 ■ 

We may apply this very simply to the case of the B(n) graphs in Figure 1. For 
each order n, the connected graphs clearly consist of a single graph. Therefore for 
each n we have c{n) — 1; whence, C{x) — exp(x) — 1. It follows that the generating 
function for all the graphs A{x) is given by 

A{x) = exp(exp(a;) - l) , (10) 

which is therefore the generating function for the Bell numbers. 

Given a general generating function C{x) — X^^i V„x"/n\ we may represent 
this by the same sequence of graphs as in Figure 1, but now with vertex factors 
Vn = c{n), as shown in Figure 1. 

This is of course a trivial observation, since any sequence of graphs each term of 
which consists of a single connected graph will, with the appropriate multipliers Vn, 
give by definition the exponential generating function of many combinatorial se- 
quences. This however shows the generic nature of the B{n) sequence, as previously 
noted in [S]; and will prove of value in our discussion of partition functions. 
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3 Partition functions 

We first illustrate the use of normal ordering in the context of partition functions 
by the most simple example. We define as usual the partition function Z associated 
with a Hamiltonian H by Z = Trexp(— 

3.1 Free boson gas and Bell polynomials 

We take as our example the Hamiltonian for the single-mode free boson gas H — 
ea^a (ignoring an additive constant), e > 0. The usual computation of the partition 
function, exploiting the completeness property X^S^Lo — immediate: 

Z = Trexp(-/3eata) = (11) 

OO 

^Y.{n\e-f'"^\n)^ (12) 

n=0 



= Y.e-^^-= (13) 

^(l~C-^'y\ (14) 



However, we may, somewhat perversely in this simple case, use any complete set 
to perform the trace. Choosing coherent states, for which the completeness or 
resolution of unity property is 



d^z\z){z\= I = J d^{z)\z){z\ (15) 
the appropriate trace calculation is 

^ = i y" d^z{z\ exp(-/3eata) \z) = (16) 

d^z{z\ : exp(ata(e-/3' - 1)) : \z) , (17) 



1 

TT 

where we have used the following well-known relation [HI d for the forgetful nor- 
mal ordering operator : f{a,a^) : which means "normally order the creation and 
annihilation operators in / forgetting the commutation relation [a, a^] — 1"^): 

Af{exp{xa^a)) exp(a^a(e^ " 1)) : • (18) 

We therefore obtain, integrating over the angle variable 9 and the radial variable 

r = |z|. 



1 p27T nOO 

Z = - dO rdrexp {r^{e-"'^' - l)) , (19) 
-'o JO 



\ Of course, this procedure may alter the value of the operator to which it is applied. 
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which gives us Z = (l — c '^^) ^ as before. 

We rewrite the above equation to show the connection with our previously- 
defined combinatorial numbers. Writing y = and x — —l3e, Ea. H19() becomes 

oo 

Z = I dyc^p{y{e- -I)). (20) 



This is an integral over the classical generating function for the Bell polynomials 
exp(y(e--l)) =5]i3„(2/)^, (21) 

where B„(l) = B{n), as defined in Eq.Q. This leads to the combinatorial form 
for the partition function 

Z= / dy^i?„(y) — . (22) 

Although Eg. 1)22(1 is remarkably simple in form, it is often by no means a 
straightforward matter to evaluate the analogous integral for other than the free 
boson system considered here. Further, it is also clear that we may not interchange 
the integral and the summation, as each individual y integral diverges. We shall 
therefore concentrate in what follows on the partition function integrand (PFI) 
F{z) = (z| exp{—PH)\z) , whence Z = J F(z) d/i(z), to give a graphical description 
of a perturbation approach. The function F maps coherent states \z) to (real) 
numbers, and is related to, but not identical with, the previously-introduced free 
energy functional ;8] associated with the Hamiltonian system. 



3.2 General partition functions 

We now apply this graphical approach to the general partition function in second 
quantized form. With the usual definition for the partition function 

Z = Trexp(-/3i7) . (23) 

In general the Hamiltonian is given by H — ew{a,a''), where e is the energy scale, 
and w is a string (= sum of products of positive powers) of boson creation and 
annihilation operators. The partition function integrand F for which we seek to 
give a graphical expansion, is 

Z{x) = f F{x,z)Aii{z), (24) 
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where 

F{x, z) = {z\ eyi'p(xw)\z) — — 



Y.{z\w-\z)- 

n=0 



X" 



n=0 



X" 



expl^KW^), (25) 



with obvious definitions of Wn and Vn- The sequences {W„} and {14} may each 
be recursively obtained from the other |Uj . This relates the sequence of multipliers 
{Vn} of Figure 1 to the Hamiltonian of Ea. (|23|l . The lower limit 1 in the Vn 
summation is a consequence of the normalization of the coherent state \z). 
We conclude this discussion with an example. 

3.3 Superfluid boson system 

We may write the Hamiltonian for a superfluid boson system as |l()j 

H^e{\ cia^ + \ ci(at)2 + C2(a^a + i) } . (26) 

As in the non-interacting boson gas example, this represents a single mode of the 
system; the complete Hamiltonian is a direct sum of non-interacting terms. The 
constants c\ and C2 are complex and real respectively, but otherwise arbitrary. 
Noting that the Hamiltonian Eq. H26|l is an element of the non-compact Lie algebra 
su(l,l) enables an expression for the partition function integrand to be readily 
obtained (see the Appendix). To illustrate the graphical series for this PFI, we 
shall choose real arguments; that is we put z and the constant c\ real in Eq.|j22I) 
below to get 

~ " 1 



F{x, y) = y i exp 



2/ yi H 1 



(27) 



where y — z^ and j/i, /i are functions of x. The PFI F{x, y) of Ea. (|27|l generates an 
integer sequence for even ci, C2 as may be easily verified using algebraic software. 
For example, choosing ci = 2, C2 = 4 we have 

Fix,y)^e.jYlMl + ^^^ + ..), (28) 



with 



1! 2! 3! 



Vi{y)=2 + 6y, 
V2{y)^2 + 36y, 

Vsiy) = 16 + 288y , (29) 



Viiy) = 144 + 3024y, 
leading to a corresponding Bell-type graphical expansion. 
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4 Discussion 

The normal ordering of creation and annihilation operators leads naturally to 
Stirling and Bell numbers. Therefore these combinatorial numbers, and their graph- 
ical representations, are ubiquitous in second-quantized quantum physics. In this 
talk we emphasized this fact by application to a perturbative expansion for the 
partition function of statistical physics. We did this through the introduction of 
the partition Junction integrand or PFI, a function which when integrated over the 
complete set of coherent states gives the partition function proper. The evaluation 
of the PFI depends on the normal ordering of the the exponentiated Hamilto- 
nian given in terms of second-quantized operators. In general this is a difficult 
if not intractable problem, although approaches based on the Product Theorem of 
combinatorial graph theory lead to some simplification'^). However, even with the 
straightforward methods illustrated here, solutions are obtainable in many cases of 
physical interest. In particular, we illustrated this approach in the case of the free 
boson gas, and a superfluid boson model. We also gave a graphical representation 
for the perturbation expansion of the PFI, an expansion in terms of a; = — /3e, es- 
sentially energy scale times inverse temperature. The graphical series is that of the 
Bell numbers, together with appropriate multipliers where z is the coherent 

state parameter which must be integrated over to obtain the partition function. 
We showed how to determine these multipliers in terms of the Hamiltonian. We 
conclude that there is a sense in which the Bell graph series is generic for the 
evaluation of PFI's. 

The operators of this theory are quantized, but have no explicit dependence 
on space or time, which justifies considering the theory as a zero-dimensional field 
theory, with the graphs as being analogous to zero-dimensional Feynman diagrams. 

Appendix: Partition function integrand for a superfluid boson system 

The partition function integrand F{x, z) corresponding to the hamiltonian Ea. H26|l 
is the expectation of the following operator 

ex.p{—PH) = exp { i xcio^ + ^ xci(a^)^ + a;c2(a^a + 5)} = 

for ^1=^^!'^ 

X2 = XC2 , y 

= exp{a;ii^- + xiK+ + 2x2K°} , (30) 

where we have identified the generators {K^, K~, K'^} of su{l, 1) with the appro- 
priate operators {^{a^)'^, ^a^, ^{a'^a + i)}. The normally ordered form of Ea. l|Sn|) 
is 

exp(-/?i/) = e^i^^e^^^^^e^^^" . (31) 

^) See the talk by A. Horzela et al in this conference. 
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Since the equality of Eq.ljSOJ and Eq.(P3 is a group-theoretical result, we may use 
the 2x2 representation 

,32) 

in which to implement the equality, thus evaluating yi and 1/2 ■ The result is 

2/1 = xci sinh ^-^^ , y2^-lnfJ,, 
where ^ ^ 

S — Xy^ C2^ — |ci p , H = cosh 6 — XC2 - 



5 

We therefore have for the Partition Function Integrand 

F{x,z) = (z|exp{ia;io^ + iii(a^)^ +2:2(0^0+ ^)}\z) = 
= <^z|e^i(''^^'/^e^"^°^°+^/^^e^i'''/^|z) = 
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